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Abstract 

The transition states and dividing surfaces used to find rate constants for bimolecular 
reactions are shown to undergo Morse bifurcations, in which they change diffeomorphism 
class, and to exist for a large range of energies, not just immediately above the critical en¬ 
ergy for first connection between reactants and products. Specifically, we consider capture 
between two molecules and the associated transition states for the case of non-zero angu¬ 
lar momentum and general attitudes. The capture between an atom and a diatom, and 
then a general molecule are presented, providing concrete examples of Morse bifurcations of 
transition states and dividing surfaces. 
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1 Introduction 


One way of finding rates of reaction is to consider rates of transport in a low dimensional 
Hamiltonian system representing the specific reaction. Some of the first examples studied using 
transition state theory consisted of bimolecular reaction in gaseous phase, i.e. 

A + B —> products, 

for two (polyatomic) molecules A and B. Provided the Born-Oppenheimer approximation holds, 
we can pass from the quantum mechanical system to a classical one, namely the Hamiltonian 
system for the motion of the nuclei interacting via a potential given by the (ground state) 
energy of the electrons 1 as a function of the internuclear coordinates. Then, if this extremely 
high dimensional Hamiltonian system is, at any instant, the product of “reacting” two molecule 
sub-systems that are independent of each other, we may consider the evolution of an ensemble 
of individual reactions in this low dimensional Hamiltonian sub-system. For this, we require the 
gas to be sufficiently dilute. See Keck [Kec67] for a nice review of transition state theory and 
discussion of the assumptions involved. Finally, since the low dimensional system is Hamiltonian 
and energy is conserved, we can restrict our attention to the energy levels, and consider the flow 
of ergode 2 , as a function of the energy. Finding (microcanonical) reaction rates translates to 
finding the rate of transport of ergode between regions representing reactants and products. 

Transition state theory provides upper bounds on rates of transport via the flux of ergode 
through a dividing surface that separates the regions of interest. The dividing surface approach, 
which we reviewed in [MS14], can be traced back to the works of Marcelin [Mar 15, Chapter 2] 
and Wigner [Wig37]. In order to obtain a useful upper bound, the dividing surface is chosen 
to have locally minimum flux in a given direction. For this to be the case, the dividing surface 
must be the union of surfaces of unidirectional stationary flux with no local recrossings, which in 
turn is the case for surfaces spanning a closed, invariant, orientable, codimension-2 submanifold 
of the energy level, known as a transition state. The simplest transition states and dividing 
surfaces are found in the basic transport scenario of “flux over a saddle” [Mac90]. These are 
transition states diffeomorphic to § 2m ~ 3 for energies just above the index-1 critical point of the 
Hamiltonian for m degree of freedom systems, with dividing surfaces diffeomorphic to S 2m-2 . 
For two degree of freedom systems, these transition states are hyperbolic periodic orbits. As 
the energy is increased, they may lose normal hyperbolicity and bifurcate. The bifurcations of 
periodic orbits are well known, see e.g. [AM78, Section 8.6], and those of periodic orbit transition 
states have been studied for a long time, see e.g. [DVB55, PP78, Dav87]. Instead, for higher 
degrees of freedom than two the transition states for the basic scenario are normally hyperbolic 
2 m — 3 spheres (within the energy levels) and very little is known about their bifurcations. 

What had been overlooked, and was explained in [MS14], is that the transition states may lose 
normal hyperbolicity by becoming singular, i.e. the manifold structure fails, at a higher critical 
energy value but then regain their normal hyperbolicity for values above the critical energy. 
These are Morse bifurcations [MS14, Appendix B], leading to a change in diffeomorphism class 
of the transition states and dividing surfaces spanning them. They occur when the union of 
transition states over a range of energies forms a smooth normally hyperbolic submanifold of 
state space, which we call a transition manifold, and there is a critical value of the Hamiltonian 
function restricted to the transition manifold. By definition, the transition states are the level 
sets of this restricted Hamiltonian and undergo a bifurcation. Similarly, the dividing surfaces, 
which span the transition manifold, undergo a Morse bifurcation. Actually, they occur because 
the energy levels undergo a Morse bifurcation themselves so the dividing surfaces, and therefore 
the transition states, must also undergo a change of diffeomorphism class in order to still separate 

*1 Assumed non-degenerate and hence a smooth energy function, without conical singularities, see e.g. [DYK04]. 

*2 Ergode is Boltzmann’s name for a microcanonical ensemble, see [Bru76, pages 242, 367]. We shall use it 

interchangeably with energy-surface volume. 
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these in two. Morse theory tells us that we can continue the transition states and dividing 
surfaces, with respect to the energy, through critical values, as well as giving the diffeomorphism 
class of the submanifolds. This is useful, for example, when the transition states above some 
critical energy are disjoint, which we shall see in the atom-frozen diatom system, as we then 
know exactly what objects to look for. Recently, there has been a lot of interest in general 
transport problems (as opposed to the basic one), non-minimum energy paths [Osb08], and 
roaming reactions [BS11] . These transport scenarios may emerge from the basic one via a 
Morse bifurcation, in which case by considering these bifurcations, we obtain transition states 
and dividing surfaces that can be used to find the rates of transport. 

We shall consider transport between the region representing two distant molecules (reactants) 
and the region in which the molecules are close. The latter does not however generally constitute 
the products 3 . This is the capture transport problem associated with the necessary first step 
of the molecules getting close enough to react. The capture rate (sometimes also called collision 
rate ) provides an upper bound on the reaction rate, as we do not expect all captured trajectories 
to proceed to the products region [CSB80]. Note that there might actually be multiple product 
regions, but for a two-body capture process there is only one final region of interest. 

Capture rates are crucial for many physical processes, and have a long history dating back at 
least to 1905 with Langevin’s early contribution [Lan05], see review by Chesnavich and Bowers 
[CB82]. Two assumptions are usually found in the literature. Firstly, the reacting Hamiltonian 
systems are assumed to have Euclidean symmetry, that is to be invariant under translations and 
rotations. This is the case for gas phase reactions with no background (electro-magnetic) field. 
The Hamiltonian system can then be reduced to a family of systems, in centre of mass frame, 
parametrised by the angular momentum. The Hamiltonian function then contains both Coriolis 
and centrifugal terms. Secondly, the energy is taken to be below those at which the two molecules 
dissociate and centrifugal and Coriolis forces to be sufficiently weak such that the molecules 
are well defined and in the small vibrations regime* * 1 . These assumptions allow to distinguish 
between intermolecular degrees of freedom (distance and relative attitudes of the molecules) and 
intramolecular ones. We too shall consider systems that satisfy these assumptions. 

We want to find the rate of capture, which we shall assume can be thought of as trans¬ 
port between regions on either side of a non-degenerate maximum x c of the effective potential 
(centrifugal terms plus potential) with respect to the intermolecular distance x. In the liter¬ 
ature, this maximum is generally a centrifugal maximum obtained by balancing the repulsive 
centrifugal terms with the attractive long distance potential energy. Alternatively, x c could be a 
non-degenerate chemical maximum of the bimolecular potential and therefore of the effective po¬ 
tential for small angular momentum. Provided x c is sufficiently large, such that capture occurs 
in a region where the potential is only weakly dependent upon the attitudes of the molecules, 
and sufficiently non-degenerate, we shall see that fixing the intermolecular distance degree of 
freedom to the maximum value gives a normally hyperbolic transition manifold in state space, 
which can be spanned by a dividing manifold. The restriction of these manifolds to the energy 
levels gives dividing surfaces and transition states, which we shall refer to as capture transition 
states. The literature often refers to them as orbiting or loose transition states [CB82, Pec76]. 

The central field model, in which the attitudes of the colliding pair are ignored, is attributed 
to Langevin [Lan05]. In this very early work, one already finds capture periodic orbit transition 
states. Langevin considers the capture process using scattering theory. For an introduction to 
scattering theory see e.g. [GPS02, Section 3.10], whereas for a comparison with the dividing 
surface approach see [CB82], Non-central fields were considered later, also usually from a scat¬ 
tering theory perspective, starting with the works of Pechukas [Pec80] and Chesnavich, Su and 

*3 Association and recombination reactions are largely limited to reactions in condensed phase or solvent, see 
e.g. [HH08, Chapter 1] 

*4 We are implicitly assuming that the molecules are normal, i.e. that their potential has non-degenerate critical 
points corresponding to rigid equilibrium configurations. Molecules that are not normal are anomalous. We 
avoid the term rigid, as it might lead to confusion with the rigid body limit with no vibrations. 
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Bowers [CSB80]. As the energy is increased above the critical energy for first connection be¬ 
tween reactants and “products”, the intermolecular attitude and angular momentum degrees of 
freedom will generally be involved in interesting sequences of Morse bifurcations causing changes 
in the diffeomorphism class of the energy levels, as the “cone of acceptance” (i.e. orientations 
for which capture is possible) opens up, and therefore changes of the capture transition states 
and dividing surfaces. Instead, the intramolecular degrees of freedom are assumed to consist of 
small vibrations, thus not playing a role in the Morse bifurcations. Using Morse theory, one can 
know which the correct transition states and dividing surfaces are, as we shall outline for two 
examples, and then study rates of transport for larger energies than previously thought possi¬ 
ble. Here, we shall focus on the Morse bifurcations, the reduction and the details of bimolecular 
capture problems. Some analysis of “reaction dynamics” of rotating molecules has been done 
recently in [CW12, KKlla], but we are interested in the interaction of two rotating molecules. 

Whether captured pairs then go on to react can be thought of as a further transport problem 
with its own ( reaction ) transition states and dividing surfaces, possibly (but not necessarily) 
associated to other maxima Xi of the effective potential. These are often referred to as tight 
transition states. The capture and reaction transition states are therefore in series. The simplest 
case will be when these are distant and the level sets of separate transition manifolds. However, 
even when “separate” their stable and unstable manifolds, which act as transport barriers, will 
intersect, determining the “reaction channels”. Some trajectories joining reactants and products 
might roam in the region between the (capture and reaction) dividing surfaces, that is follow 
trajectories with a non-monotonic intermolecular distance in time, before finally crossing the 
reaction dividing surface. This is to be expected due to coupling between degrees of freedom, 
and was recently given as an explanation of what has been called roaming reactions [MCE+14], 

Reaction rates have an equally long history as capture rates, and bimolecular reactions 
have played the role of test problems since the early days of transition state theory (as noted 
in [Wig38]). The transport problems associated with reaction tend to be harder because the 
chemical potentials are at best not simple and often degenerate. Similarly to how the first capture 
models were simplified by making the fields central, reaction rates were first, and largely still, 
considered for collinear and planar systems with zero angular momentum. Note that reaction 
transition states will also generally undergo Morse bifurcations. 

Hamiltonian systems with symmetry, such as the molecular n body systems with Euclidean 
symmetry considered here, have conserved quantities by Noether’s theorem (see e.g. [Arn89, Ap¬ 
pendix 5], [Mar92, Section 2.7]), and can be reduced via symplectic reduction. The Euclidean 
symmetry group of our n body systems is a semi-direct product SE(2>) = M 3 x SO(2>) and can 
be reduced in stages [MMO + 07, Chapter 4], The translational symmetry and associated linear 
momentum are easily reduced, whereas the rotational symmetry, which does not act freely on 
the whole of state space, and the associated angular momentum require singular reduction, see 
e.g. [SL91], [OR04, Section 8]. The reduced state space is therefore a stratified symplectic 
manifold, however the Hamiltonian flow leaves the connected components of the strata invariant 
so the reduced dynamics can be studied on each stratum individually. The stratified symplectic 
manifold for our reduced molecular n body systems is composed of three strata. Both singular 
strata have zero angular momentum, whereas on the principal stratum the angular momentum 
is arbitrary. This is discussed in [Strl5]. Actually, reactions with initial positions and momenta 
confined to a plane remain in this plane for all successive times, and reduced planar systems 
have no angular momentum degree of freedom and no coordinate singularities at collinear con¬ 
figurations. In Section 2, we therefore consider planar capture between an atom and a diatom 
and find our first examples of Morse bifurcations. This example and that of two interacting 
diatoms were considered in less detail in [MS 14]. A spatial example, of an atom interacting with 
a molecule, is considered in Section 3. We choose to use the coordinates that one obtains via the 
gauge theoretic approach to cotangent bundle reduction as outlined by Littlejohn and Reinsch 
[LR97], and reviewed in Appendix A, where we discuss the symplectic form and consider the 
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Figure 1: Choice of Jacobi vectors and reduced coordinates q = (x,b, /3) for planar atom-diatom 
reactions. 

transformation to Serret-Andoyer coordinates on the angular momentum sphere, which is not 
done in [LR97]. This approach gives coordinates, which are furthermore physically meaningful 
since coordinates and momenta are not mixed, however it involves considering the rotational 
50(3) action on configuration space as opposed to the lifted action in state space, and therefore 
introduces coordinate singularities. These are avoided by considering non-collinear molecules. 
Ways to deal with them are discussed in the Conclusions in Section 4. 


2 Planar atom-diatom reactions 


The simplest non-trivial capture example is planar reactions between an atom and a diatom. The 
planar reduced three-body Hamiltonian system that represents an atom interacting with a di¬ 
atom with no background field is the family of mechanical systems 5 oj, H ), parametrised 

by the angular momentum A € K, with 


H(z-X) = 


1 2 1 2 

—Px H Pb + I ,2 

m mb \mi>b z 


1 1 

+ 


mx- 


P0 


rribb 2 


mx 2 + 7Tibb 2 


A + V(q;\), 


A 2 


2 mx 2 + myb 2 


+ U(q) 


where z = (q,p), q = (x, b,j3), x is the intermolecular distance, /3 the attitude and b the length of 
the diatom, see Figure 1. V is the effective potential with the centrifugal term. The parameters 
are the reduced masses m and nib, and the magnitude of the angular momentum A. 

We shall restrict our attention to energies below that at which the diatom dissociates, so we 
have a two-body capture problem. Specifically, we shall consider the capture of an atom by a 
strongly bonded diatom that is distant and rotating slowly, assuming that the potential of the 
pair has a non-degenerate maximum x c with respect to x for large x, and that the centrifugal and 
Coriolis forces are not too large. Thus, for energies below that at which the diatom dissociates, 
the reduced coordinates and their momenta split into intermolecular (x, /3) and intramolecular 
b degrees of freedom. Our assumptions imply that the diatom is in the small vibrations regime, 
i.e. ( b,pb ) is an elliptic degree of freedom, and rotating slowly, so in the neighbourhood of x c 
we will find a codimension-2 normally hyperbolic submanifolds N and span it by a dividing 
manifold 5. This will be the capture dividing manifold that trajectories must cross in order for 
capture to occur. The assumptions will be introduced by scaling the coordinates, after which 
the dynamics and the existence of the transition manifold N will be clear. Then, we will show 
that the level sets of N and 5, the transition states Ne and dividing surfaces Se, generically 
undergo Morse bifurcations involving the (/3,pp) degree of freedom, as the energy is varied. 

For simplicity, the masses of the atom and the diatom are set to m = mb = 1. A better 
approach would be to non-dimensionalise the variables. Note that normalised or mass-weighted 
Jacobi coordinates do not remove the mass dependence. 


*5 A simple mechanical system is a Hamiltonian system whose state space is the cotangent bundle of a Rieman- 
nian manifold ( configuration space) with canonical symplectic form, and the Hamiltonian is the sum of the 
positive definite kinetic energy, given by the metric, and a potential energy. 
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Figure 2: Left: Typical graph of molecular potential restricted to the intermolecular distance, 
with repulsive short range, attractive long range and extrema in between. Right: Example 
graph of the effective potential over the intermolecular distance and angular momentum (x, A) 
showing the disappearance of the intermolecular maximum x c , via centre-saddle bifurcations. 


The intermolecular terms of the potential will be repulsive at short ranges, possibly have a 
number of maxima (and therefore minima) with respect to the intermolecular distance x in the 
mid ranges, and be attractive at long ranges with inverse fc-power of x terms [Stol3], see Figure 
2. The potential is then summed to the repulsive centrifugal term to give the effective potential. 
In the long (physical) range, provided the attractive potential falls off faster than the centrifugal 
potential as x —> oo, i.e. k > 2, the effective potential has a centrifugal maximum x\(b, /3; A). In 
the short (chemical) range of the potential, the chemical maxima of U with respect to x imply 
chemical maxima of V for A small with respect to the slope of U at the maxima. In either case, 
as A increases the maxima will “collide” with the minima and disappear, see Figure 2. 

We shall think of capture as transport between the regions on either side of the largest 
maximum x c (b,/3]\), assumed large with respect to the length of the diatom b. We introduce 
the capture scale by setting b = e c b and p b = e^Pb with 0 < e c « 1, and taking x c to be of 
order 1. Essentially, e c is the ratio between the size of the diatom and its distance to the atom, 
however in practice we assume that for distant pairs, the potential energy is weakly dependent 
on the attitude of the diatom, since the pair are distant, and we can choose e c such that 

U(q-e c ) = e~ 2 U b (b) + U°(x) + s 2 c U 2 (q- e c ). 


This is the case for potentials that are inverse power functions of the intermolecular distance, 
which can be expanded using Legendre polynomials. The Hamiltonian function expands into 


2 

-c 


jrf + ifl + EW)' 


H(z; A, e c ) = £ 

+ e lU 2 (q] 0) + O (ec) > 


+ 2 W 


1 2 , 1(P0~A)" 


+ 


+ U°(x) 


where we have dropped the tildes. We note that f3 does not appear until order e^, so pp = 
A p + O (e^) with constant A^. Furthermore, the system separates into slow and fast degrees of 
freedom, i.e. at order e ~ 2 we find the fast oscillations of the diatom plus a “centrifugal” term, 
at order there is the intermolecular (capture) dynamics, and then there are the higher order 
terms. Up to order the x and b degrees of freedom are uncoupled, and pp = Xp. 

The diatom will have an equilibrium configuration, which corresponds to a non-degenerate 
minimum b(x, /3) of the joint atom-diatom potential with respect to the intramolecular distance 
b. We assume this minimum to be highly non-degenerate, i.e. that the diatom is strongly bonded, 
and that the centrifugal and Coriolis forces on the diatom are not too strong, so the diatom will 
vibrate about its equilibrium without significant distortion. We therefore, linearise b about 


6 







by setting b = bo + Sbb and pb = £ b l Pb with 0 < £& <C 1. The constant £b is chosen such that 

U{q ;£) = £ c 2 (&b + 2 £ b 2 Ubb 2 'j + U®(x) + £ 2 U 2 (q;£ c ) + O (e^) , 

where U b := t4(&o), U‘b := £ b d 2 h Ub(bo) is order one, and £ = (e c , £■&). Recall that we are assuming 
d‘hiPb(bo) to be large. This scaling ensures that the leading order terms of the potential with 
respect to the coordinates are of the same order as their conjugate momenta in the kinetic 
energy. The Hamiltonian function, again dropping the tildes, becomes 

H{z- X, £ ) = £c 2£ b 2 \ (.Pb + Ufa 2 ) + £ c 2 ^2 p} + (\pI + 7^ (P8 ~ + U c( X )^j 

+ £ c^ 2 (bo,x, /?; 0) + O (^c, £ l) ■ 

We shall further simplify our Hamiltonian by setting pp = £^.pp, i.e. considering \p = 0. 
General pp is considered in the disconnecting example of [MS 14]. The scaled system consists of 

H(z\ A, e) = £ c 2 £ b 2 - (pi + U 2 b 2 ) + -p 2 + + UcW 

+ e c + ~^PP X + Uc(b o,x, (3-, 0)) + O ( £c, £ l) 

and 

lo = db A d pb + £ 2 df3 A dpp + dx A d p x , 
from which we get the equations of motion up to order e°, namely 


b £c £ b ~Pbi 


Pb £c £ b 


1 1 

P = 79 P0 - 2 A ’ X = P*' 

°5 x 

P/3 = -dpU 2 (b 0 ,P,x-0), p x = ^3 A 2 -d x U°(x). 


By assumption, the intermolecular distance degree of freedom is hyperbolic, and the in¬ 
tramolecular distance is in the small vibrations regime, i.e. elliptic. These dynamics are uncou¬ 
pled to this order. As the diatom rotates, the attitude degree of freedom will display both kinds 
of motion. Provided the (x,p x ) degree of freedom is more strongly hyperbolic that the ( i 0,pp ) 
one, that is the maximum x c is sufficiently non-degenerate, the submanifold 

No = {z G M\\x = x°(A), p x = 0} 

is almost invariant, or more precisely invariant to order e°, and normally hyperbolic. By normal 
hyperbolicity theory, there is a true normally hyperbolic submanifold N nearby. Given No, we 
could find a better approximation to N as explained in [MS14], However, for the purpose of 
studying the Morse bifurcations of the transition states, No is a sufficiently good approximation. 

The normally hyperbolic submanifold N is a transition manifold, as it can be spanned by a 
dividing manifold. The approximate transition manifold Nq is spanned by 


So = {z G M\\x = x c (b, /3; A)}. 


The transition states Ne and dividing surfaces Se are then approximately the level sets of the 
Hamiltonian restricted to the approximate transition and dividing manifolds. As the energy 
varies, we expect these to bifurcate. The transition states may lose normal hyperbolicity. For 
atom-diatom reactions, Ne are 3-dimensional manifolds, so it is not well understood how they 
lose normal hyperbolicity. Freezing the diatom, the system only has two degrees of freedom and 
Ne is a periodic orbit. In [MS14], we saw that these disappear in a centre-saddle bifurcation. 
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Figure 3: Graphs of an example frozen Hamiltonian restricted to the transition manifold H^ 
(right) and the potential U 2 (left) for an atom-diatom reactions with the atom attracted to one 
of the sides of the diatom (e.g. ion plus dipole), with intermolecular and diatom distances fixed. 


However, before the loss of normal hyperbolicity, the capture transition states will undergo 
changes of diffeomorphism class via Morse bifurcations. 

If we write 


N = {z € M\\x = 4(A) + 44(4 + O (4) , Px = 0 + 4^M4 + O (4)}, 


we find that the Hamiltonian function restricted to the transition manifold N is independent of 
x 2 and P 2 up to order 4, namely 


H n (z;X,s) 



(p 2 b + U 2 b 2 )+s 2 c 



+ E 2 (b 0 , Xq, 



52 

where vr = pn — -QX + • • • is the non-canonical momentum, and we dropped constant terms. 

H r x o 

For E below that at which the diatom dissociates, the intramolecular degree of freedom 
contributes only positive-definite terms to the restricted Hamiltonian function and is not involved 
in any Morse bifurcations, which can be studied by considering the simpler function obtained by 
freezing the diatom, i.e. minimising Hjy over (b,p b ) by setting b = 6q + h.o.t. and p b = 0 giving 


H$r(P,pp\X,e) = 4 



Different reactions, i.e. choices of atom and diatom, will have different potentials and different 
sequences of Morse bifurcations. An example frozen restricted Hamiltonian is depicted in 
Figure 3. Considering this H^, we find that N< E = (Ff^-) _1 ((— 00 , E}) bifurcates from B 2 to 

S 1 x B 1 , and so the level sets N E from S 1 to §° x S 1 , where B denotes a ball and § a sphere. 
Passing to the full system, we find that the transition manifold N<e bifurcates from B 4 to 
S 1 x B 3 , the transition state Ne from § 3 to S 1 x § 2 , and the dividing surface Se from § 4 to 
S 1 x § 3 . Thus, the energy levels are separated by the dividing surfaces for all energies, even 
above the critical one of the Morse bifurcation. Other examples of U 2 are considered in [MS 14]. 

Given a dividing surface Se, we can compute the flux of ergode through it, denoted 4>e{Se), 
and use this to find the rate constant. The flux varies C m ~ 2 smoothly through the Morse 
bifurcations, for m degrees of freedom, but not C m ~~ l [MS14]. Graphs of cj>E as a function of E 
for the capture dividing surface of the atom-diatom system with the example potential presented 
in Figure 3 can be seen in Figure 4. For the periodic orbit transition state of the two degree 
of freedom system with the frozen diatom, we see the log-like infinite slope singularity at the 
index-1 Morse bifurcation. This is not present in the graph of the flux through the dividing 
surface of the three-degree of freedom system. 
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E E 

Figure 4: Example graphs of the flux 4>e as a function of energy E, for the capture transition 
state of the atom-diatom system, with potential as in Figure 3, undergoing a disconnecting 
Morse bifurcation at E 2 . Left: two degree of freedom system with frozen diatom. Right: three 
degree of freedom system. 


3 Spatial atom-molecule reactions 


We now consider the capture of an atom by a polyatomic molecule. The molecule shall be 
assumed to have a non-degenerate non-collinear equilibrium, of its rib atoms, about which it 
is vibrating fast, and the system to have energy below that at which the molecule dissociates, 
so we can use the charts provided by the bundle approach to cotangent bundle reduction, see 
Appendix A. The capture scenario will be the same as the atom-diatom one, modulo differences 
due to its spatial nature, such as an additional angular momentum degree of freedom. We shall 
again assume that the effective potential has a non-degenerate maximum x c with respect to the 
intermolecular distance x when the pair is distant, that the attitude terms are not dominant in 
this range, and that B is a normal molecule with a very non-degenerate equilibrium, i.e. that it is 
strongly bounded. In order to have B in the (elliptic) small vibrations regime, and the attitude 
and angular momentum degrees of freedom dominated by the hyperbolic (in the neighbourhood 
of x c ) intermolecular one, we ask that the molecule is rotating slowly with respect to the atom, 
that it is 3D with distinct principal moments of inertia and that the energy is bounded such 
that, roughly speaking, the angular momentum l is almost perpendicular to the line of centres. 
These extra constraints due to the spatial nature shall be discussed in detail after we scale the 
intermolecular and intramolecular distance. Finally after also scaling the angular momentum 
coordinates and the attitude momenta, we shall find a normally hyperbolic transition manifold 
N and dividing manifold S about x c and discuss the possible Morse bifurcations of the transition 
states Ne and dividing surfaces Se- These will now also involve the angular momentum degree 
of freedom, degeneracies when the moment of inertia tensor has equal moments, and require an 
understanding of the diffeomorphism class of the reduced space, which we shall discuss. 

Capture between two non-collinear molecules can be studied following the same steps as 
those below, and will be presented in a future publication. 

We denote the reduced molecular n = rib + 1 body Hamiltonian system by H\). 

The state space M\ is the subset of the principal stratum with non-collinear configurations. It 
is a smooth (6 rib — 4)-manifold diffeomorphic to a (generally non-trivial) § 2 fibre-bundle with 
base space T*(Qi d /SO( 3)), where Qj d is the non-collinear subset of the translation reduced 
configuration space Q = R 3 ( n_1 ) ; and fibre the angular momentum sphere S| [MMO + 07, Section 
2.3], 
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In canonical coordinates z = (q, q\,p,p\) £ M x , we have 

.. 3n—6 3 

H{z;\) = - Y Y^ Pi - A ik{q)h{z\\ (q)(pj - A jk (q)l k (zy, A)) + V(q,zy, A) 

i,j=1 k= 1 

1 3 

V{q,Z\;\) = - Y li(z x ;\)r\q)l j (z x -,\) + U(q). 
ij— 1 

where q are the internal coordinates, which we shall choose shortly, A the magnitude of the 
angular momentum and z x = ( q X ,P \) canonical Serret-Andoyer coordinates on the angular 
momentum sphere, such that e.g. the angular momentum is given by 

= {p\, \J X 2 — p\ sin q X ,\j X 2 -p 2 x cosq x ). 

More than one chart is necessary, due to the coordinate singularities both for Qj d /SO( 3), with 
n > 4 bodies, as well as on the angular momentum sphere. V is the effective potential with 
the centrifugal terms, K is the reduced metric, I is the moment of inertia tensor, and A is the 
gauge potential, present in the Coriolis terms. These are introduced in Apprendix A. We are 
uninterested in scaling effects due to the mass, so we set mi = 1. This affects I\ and /. 

The rotating frame for the reduction is chosen such that the Jacobi vector along the line of 
centres r nb (q) is parallel to the x\ axis, and the remaining 50(2) symmetry about the aq-axis is 
used to orient the equilibrium configuration of the molecule B such that its moment of inertia 
tensor I b = Diag(/Ltf,i, p b 2 -,Pb 3 ) with p b i > p b2 > p b 3 . This determines the reduction of the 
50(2) symmetry about x\ and the definition of the attitude of the molecule. The most natural 
choice of reduced coordinates q is the distance between atom and molecule x and two angles 
[3 = (/3i,/3 2 ) £ 50(3)/50(2) = § 2 for the attitude of the molecule, which are intermolecular 
coordinates, as well as 3(n& — 2 ) coordinates b for the intramolecular degrees of freedom of 
B, so q = (x,j3,b), unless B has further symmetries of its own that can be reduced. The 
intramolecular coordinates b shall be chosen in order to simplify the Hamiltonian along the lines 
of the Eckart [Eck35] and Sayvetz [Say39] conventions for normal and anomalous molecules in the 
small vibration regime. Essentially, we shall consider an Eckart convention for a normal molecule 
in the small vibrations regime interacting with an atom, for which the intermolecular coordinates 
are similar to the large amplitude coordinates of anomalous molecules considered by Sayvetz (see 
Apprendix B). In the Eckart convention the rotations and vibrations are decoupled to leading 
order since the intramolecular coordinates b are chosen to be Riemann normal coordinates for 
which the gauge potential A\ } (q) vanishes at the equilibrium configuration. This is discussed 
from a geometric perspective by Littlejohn and Mitchell [LM02]. 

First, we scale the distances and their momenta. We are interested in large x about the (cap¬ 
ture) maximum x c , so we scale x = and p x = e c p x - Then, by passing to the intermolecular 
time, we scale the Hamiltonian such that x is of order one. Also, we consider the molecule in 
the small vibrations regime about its non-degenerate minimum, and shift the intramolecular 
coordinates to have b = 0 at equilibrium, and then scale b = £\,b and pi, = £' b 1 Pb■ As for the 
planar case, we assume that the potential scales to 

U(q;e) = U b (b) +e 2 c U°(x) + £ 4 c U^(q;e), 
and then choose £ b such that 

3ri(,—6 

£% e) = B° + £ b 2 Y Ufa 2 + e 2 c U° c (x) + £ A c U 2 c {q, e) + O (eg) . 

2—1 

That is, we are assuming that the molecule is strongly bonded, d% b Ub( 0) ~ £ b 4 . Note that we 
have chosen normal mode intramolecular coordinates for which The reduced 
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Figure 5: Angular momentum sphere with equipotential lines of the centrifugal energy, when 
the moment of inertia has three distinct principal moments (left), or two equal moments (right). 


kinetic and centrifugal energies contain both intermolecular and intramolecular terms and are 
scaled in Appendix B. The Hamiltonian function becomes 


--2 


3 ri £,—6 


H («; A > £ ) = -y- \ e b 2 U bi(Pbi + H A ’ £ ) J Jo (P) v Pi( z 'i £ ) + J 0 n (^)PA 


2 

2 1 3 1 2 
+ ^ 9 1+ 9 A) + U°(x) - - ^ vpi(z; \,e)J^ 2 (x, (3)vpj(z-, A,e) 


Z=1 


»j=l 


i,J=l 


M=1 


+ y | X] \z)J l ih(. x ,P)vpj{z- 1 \,E) + ^( z a; A) + U 2 (q;0) 

\*,i=i *,i=i 

+ h.o.t., 


where I\p 1 (q) = K^(/3)-e 2 c J 2 (x, f3)+ei J 4 (x, /3) and I 1 (q) = I Q 1 (/3)+s 2 c I 2 1 (x,^)+e^I 4 1 (x,/3) 
for some Jj, see Appendix B. We are using angular momenta l and the non-canonical momenta 

3 

v/3i(z-, A,e) =ppi - ^^Apijix,/3;e)lj(zy, \), 
j = 1 

where Apij(q) = (Apn(l3), £ 2 c Api 2 {x, /3),slAp i3 (x, /3)) H-, as place-holders. 

Next, we consider the angular momentum degree of freedom z\ € From the scaled 
Hamiltonian, we note that its dynamics are coupled to the internal one and that it appears in 
terms of orders both e ~ 2 and e° c . In order to ensure that we have a normally hyperbolic capture 
transition manifold about x c , we need that the intermolecular distance degree of freedom (x,p x ) 
is more hyperbolic than the angular momentum degree of freedom. We will therefore find the 
necessary scaling for this to be the case and apply it. 

Considering the atom and the molecule as constituting a single “body”, if the molecule is in 
equilibrium with itself and with respect to the atom, we obtain a rigid body and the Hamiltonian 
reduces to the centrifugal terms. Rigid bodies follow closed curves on the angular momentum 
sphere with equilibrium points when l is parallel to the eigenvectors of / -1 , called principal 
axes [Dep67]. Typical rigid body dynamics for the case of distinct eigenvalues, or principal 
moments , is depicted in Figure 5. In general, the angular momentum degree of freedom doesn’t 
follow closed curves on S|, since it is coupled with the internal “deformation” dynamics. The 
equilibria of the system occur when q is a critical point of the effective potential V, v = 0, and 
again l is parallel to the principal axes. We therefore consider the centrifugal energy 

E x = h T (z x -X)r\q-e)l(z x -X), 
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which is not conserved, in more detail. 

First, we consider the moment of inertia tensor and its principal moments and axes. The 
eigenvectors rji of I~ 1 are the same as those of I, whereas the eigenvalues m are the reciprocals. 
I has real eigenvalues, since it is a real positive-definite matrix, and if these are distinct then 
the eigenvectors are orthogonal. To order e®, the eigenvalues and eigenvectors of / _1 are 

= I/fiP) + ■ ■ ■ > m (?) =xi + 0 (e 2 c ) , 

m{q) = e\x~ 2 + E^X^-\ -, rji(q) <G {x 2 ,x 3 } + O (e 2 ) , 


where 

m(/?) = ^(4(0) + 4(0)) ± ^(4(0)-4(0)) 2 +4(0) 2 , 

and i = 2,3. Thus in order to have distinct principal moments and axes, we require either 
4(0) + 4(0), or 4 08) 7 ^ 0. For our 3D molecule B with three distinct moments, as it 
rotates relative to the distant atom, we expect to find three pairs of points (on the attitude 
sphere S^) at which two of the moments of I(q), and so I~ 1 (q), are non-distinct. 

Then, we consider the actual centrifugal energy and its critical points. Given a fixed config¬ 
uration q with distinct p\{q) > H 2 (q) > P 3 {q), we write 


l(q, z \] X) = P\Vi (<?) + \/A 2 - P 2 sin Q x r) 2 (q) + \JX 2 - P 2 cos Qxm ( q ), 

i.e. consider Serret-Andoyer coordinates Z\ obtained by projecting onto the principal axes, 
cf. Appendix A. Then 

E \ = \ + /h 2 (<?)(A 2 - P 2 ) sin 2 Q x + /x 3 (g)(A 2 - P 2 ) cos 2 Q x ) , 

so the critical points are Z\ = (0,0), (0,7r), Z\ = (§,0), (4^,0), and Z\ = (Qa,±A). The 
superscript denotes the principal axis to which l is parallel. For non-distinct eigenvalues p 2 (q) = 
H 3 (q), choosing generalised eigenvectors for 772 , 773 , gives 

E\ = ^((Mi(?) - 2/i2 {q))P x + 2 ^ 2 (g)A 2 , 

so the critical points are the degenerate (Qa,0), and (Qa,±A). This is depicted in Figure 5. 

For arbitrary A ~ 1, in order to ensure that the centrifugal terms are higher order than the 
intermolecular (. x,p x ) ones, we restrict our attention to energies just above the critical 


E l = yA 2 = \{e 2 c x 2 + V24(0)) + ■ ■ ■ ~ e cA 2 

such that the angular momentum degree of freedom is confined to a small annulus A 2 that 
doesn’t contain Z\. If all the energy of the system is in E\, then for the non-distinct case we 
have Pa = 0 at E x = E 2 , whereas for the distinct case 


P'x(^- E l) = 


\ A 2 (7724 — ^ 34 ) 

x 4 qii 


+ ... 


4 \ 2 
~ £ n A . 


Furthermore, px = Px + O (e 2 ), so bounding E < E 2 + A, with A small, gives px = 0 + O (e 2 ) + 
O (4). Thus, we scale p x = £ 2 0A, giving 

l = A(0, sin q x , cos q x ) + £ 2 px(l, 0, 0) + O (4) . 


Next, we consider the rotational degree of freedom, with non-canonical rotational momenta 
v/ 3 i(z- A, e) = ppi - e 2 (. Apn(P)px + Ap i2 (x , /3)A sin q x + Ap i3 (x , /3)A cos q x ) H- 
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and ensure that the ( x,p x ) dynamics dominates this one also. Even though we have removed 
the coupling of vibrations and rotations up to order e°, ppi are not conserved since the reduced 
metric K is a function of p. However, the rate of change of pp is a function of p 2 p up to order 
e 2 , so we consider a molecule that is initially rotating slowly, which will take a long time to 
increases its rotational velocity. Thus, we scale ppi = £ 2 ppi, so vpi ~ e 2 . 

Finally 


H (z; A, e) = e c 2 e 


3n i>— 6 f T 2 i \ 2 

2 'b 2 + 4) + 2 P * + 2x2 + 

1=1 


X 


+ 4 I \ vp i (z-X,0)J^ o (P)vp j (z-,\,0) + \ll l {P)pi 


i,j= 1 


3 3 \ 

+ t / P X I 1 2 \x,P)l^ q y,X) + 1 - l°{qx;X)ll 3 (x,/3)l°(q x ;X) + U 2 (q;0) + 0(4,4) 


3 =2 


i,j=2 


and 

3 n (,—6 2 

uj = ^ dbi A d p u + e 2 c ^2 d A A dp/3i + £ 2 c dq x A dp x + dx A dp x , 

i= 1 2=1 

giving the equations of motion, up to order £°, as 

k = £~ 2 £^ 2 UuPbi, Pi = d Pfji H 2 (z- A), q x = d Px H 2 (z ; A), x = p x , 

Pb = — £~ 2 £p 2 U'inbi , ppi = -d i6 ,H 2 (z; A), PA = -d qx H 2 (z; A), p x = -^V^x; A), 

where H 2 (z; X) ~ e 2 . Thus, about the maximum x c , the (x,p x ) degree of freedom is more 
hyperbolic than both the attitude (P,pp) and the angular momentum z x ones. The submanifold 


N 0 = {z e M x \x = x c (A), p x = 0} 

is therefore invariant to order £°, and normally hyperbolic. Taking Nq as an approximation to 
the true normally hyperbolic submanifold N nearby, and 


So = {z € M a |x = x c (A)}, 


as an approximate dividing manifold, we can find the restricted Hamiltonian functions 

3 n b 6 rr2 / -i 3 

H n {z- A, e) = e~ 2 e^ 2 J2 -y(Pbi + 4) + 4 f - v p4 z '’ A > d)G l p 0 {P)vpj{z] A, 0) 

2=1 ^ 2 , j=l 

3 3 

+ \ I l 1 iP)P\ + ^2pxI 2 (x c ,P)l < j{q x -,X) + 1^2 ^°(«a; A)/l J (x c ,/3)^°(g A ; A) 

j=2 i,j =2 

+ U 2 (q; 0)) + O (4>4) > 

modulo constant terms, and Hs to leading orders. These give the transition states Ne = Hjj l (E) 
and dividing surfaces Se = respectively. 

As for the planar examples, it is simpler to study the Morse bifurcations if we minimise the 
reduced Hamiltonians over the positive-definite coordinates, namely b,pb, vp and p x . This can 
be simplified by using canonical angular momentum coordinates Z x aligned with the principal 
axes, and setting b = pb = vp = P x = 0 in Hjy, to obtain 

Vp (P, Q\; A, £) := H n (x c (A), 0, P, Q x , 0,0,0,0; A, e) 

= 4 (y (+2i(P) sin 2 Q x + P3i(P) cos 2 Qx) + t7 2 (/3;0)^ +0(4,4) ■ 
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We are then interested in the level-sets of and their bifurcations. The domain of is a 
subset of N, which is codimension-2 in the reduced state space A 1\. The latter is some fibre 
bundle over T*(Qj d /SO( 3)), so also N and S could in general be non-trivial bundles. We shall 
restrict our attention to subsets for which the bundle is trivial. Furthermore, we are considering 
energies below that at which the molecule dissociates, and up to just above the centrifugal energy 
for the angular momentum aligned with the 772 (/3) principal axis. 

The critical points ( (3 , Qx) of the frozen, restricted effective potential Vfc are given by 

{H 2 a(P) - fJ-te(P)) sin Qx cos Q\ = 0, 
y {dpmW) sin 2 Q\ + dpfj, :i4 (f3) cos 2 Qx) + dpU^(/3 ; 0) = 0. 

The first equation is satisfied trivially for (3 at which the two principal moments are equal 
M24(/3) = h 34 (/d). We shall consider examples of V ^ that are Morse functions, i.e. have non¬ 
degenerate critical points ( f3 , Qx) with (3^/3 , satisfying either 

Ql = kn and dp ^y/r 3 4 + U^j (/3) = 0, 

01 2 

Ql = kn + ^ and dp ^y/r 2 4 + U^j (/3) = 0, 

for k £ Z, cf. [LR97, Section IV.E], Furthermore, the Morse function has at least two 
non-degenerate minima at ( f3° , Q 3 ) due to the symmetry of the centrifugal terms. 

The sequence of Morse bifurcations of the level sets of Vfa, and therefore of Ne and Se, 
depends on the relative size of the centrifugal and the reduced potential V 2 energies. This will 
determine the relation of the different critical energies. Critical points with the same attitude 
j3 but the angular momentum aligned with different principal axes have energies that differ by 
A 2 , whereas the difference in energy for different attitudes depends on the atom-molecule pair. 

The simplest case is when the first Morse bifurcation encountered as the energy is increased 
from the minima involves the angular momentum angle, and the system goes from rotating 
about the 773(9) axis to rotating more freely about 772 ( 9 ) as well. This bifurcation occurs at 
the critical energy for the (/3°,0 a) critical points. In this case both the domain of Vfp and the 
subsets of N and S of interest are bundles over a contractible base space and so trivial [Ste51, 
Theorem 11.6]. The frozen energy levels N<e bifurcate from §° x B 3 to S 1 x B 2 , so the transition 
states Ne go from §° x to S 1 x § 6n6 ~ 8 , and similarly the dividing surfaces. 

As the energy is increased further, we will reach critical values at which also the attitude 
coordinates are involved in Morse bifurcations. We will consider the case in which the energy 
does not change significantly as the molecule rotates in one direction, with respect to the atom, 
but does when it tries to rotate in the other direction and a potential V 2 on § 2 that has a 
minimum /3°, a saddle /3 1 and two maxima (3 2 , restricting our attention to the annulus A 2 C § 2 
containing /3° and /3 1 . Therefore, choosing the angles appropriately, only one is involved in Morse 
bifurcations, whereas the other contributes positive definite terms. The subset of the transition 
manifold N of interest is a bundle over S 1 xB 6 ™ 6-9 which we claim is trivial. Firstly, we note that 
it is equivalent to the product of a bundle over S 1 with B 6 " 6-9 via homotopy-type arguments 
[Ste51, Theorem 11.4], cf. bundles over contractible spaces being trivial. The characterisation of 
bundles over spheres with structure group G depends on certain homotopy groups of G [Ste51, 
Theorem 18.5]. Our fibres are diffeomorphic to § 2 , or subsets of it, and the diffeomorphism 
group of § 2 is the orthogonal group 0(3). However, N is orientable so both elements of the 
product must be orientable. Thus, given that the bundle over the circle is orientable, we restrict 
our attention to the orientation preserving diffeomorphisms 50(3) and find that the bundle over 
the circle is a product, and therefore our original bundle is also trivial [Ste51, Section 26]. Note 
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Figure 6 : Contour plots of example functions on the (/?i,Qa) torus, where darker regions rep¬ 
resent lower energies. Top row: for no values (3\ with non-distinct principal moments. Bottom 
row: for non-distinct principal moments [12 and ^3 at depicted by vertical back lines. Case 
in which the value of the function at (j3%,Q\) is smaller than that at 011 the left, and 

vice-versa on the right. 
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however that not all orientable surface bundles over the circle are product bundles, as we can 
construct non-trivial bundles with fibres diffeomorphic to T 2 , for example. 

Vft can be minimised over the irrelevant attitude to obtain a function on the torus T 2 for 
(Pi,Q\), say. There are two possible scenarios for this case, the first is that T 2 does not contain 
points /3 at which the /J2, /U3 principal moments become equal. The order of the bifurcations 
then depends on the relative heights of the critical energies, and both cases are straightforward, 
see Figure 6. The other scenario is when T 2 does contain j3. We shall consider the case in which 
it contains only one pair of such points. Contour plots for the restricted function on T 2 are given 
in Figure 6. If the centrifugal energy is smaller than the attitude potential, then the points j3\ 
at which the moments ^ 2 , ^3 are not distinct do not play a role in the Morse bifurcations, which 
are the same as those for the case when T 2 does not contain /3, as we can see by comparing the 
left hand side of Figures 6. Instead, when the molecular potential is smaller than the centrifugal 
one, depicted on the right in Figure 6, we see that the points f3\ do play a significant role in the 
bifurcations and the sub-level sets of the torus bifurcate as follows 

§° x B 2 to §° x §° x B 2 to S 1 x B 1 to A c to Y c to T 2 , 

where X c and Y c can be written as handlebodies using Morse Theorem B [MS 14]. Therefore, the 
sub-level sets of the capture transition manifold N<e have the following sequence of bifurcations 

§° x B 6nb_6 to §° x §° x B 6Tlb_6 to S 1 x m 6nb ~ 7 to X to Y to T 2 x B 6ni >" 8 , 


and the transition states 

§° x S 6 ^" 7 to §° x §° x B , 6nb ~ 7 to S 1 x S 6,lb ~ 8 to dX to dY to T 2 x S 6 ^” 9 . 

Similarly for the dividing surfaces. 

Finally, if we were to consider higher energies for this choice of t/ 2 , the other attitude would 
also become involved in Morse bifurcations. Here again the j3 points would most likely lead to 
interesting sequences of Morse bifurcations, however we would also have to deal with the non¬ 
trivial nature of the fibre bundle. After the Morse bifurcations at the index-2 critical points /3 2 , 
the base space would contain a 2-sphere, and many examples of non-trivial orientable bundles 
over these can be found. Thus before we can consider the full sequence of Morse bifurcations of 
the dividing surfaces and transition states and the transport for a larger range of energies, the 
bundle class of the reduced state space needs to be understood. 

4 Conclusions and discussion 

The purpose of this article was to show that Morse bifurcations must be considered when study¬ 
ing transport problems for larger ranges of energies, and more specifically to show the existence 
of Morse bifurcations of capture transition states and dividing surfaces for bimolecular reactions. 
By considering the different sequences of Morse bifurcations we were able to find interesting new 
transition states and dividing surfaces for general reactions with non-zero angular momentum, 
thus extending the dividing surface approach beyond the well known basic transport scenario. 
Other choices of molecules than those considered here will have similar capture transport prob¬ 
lems and therefore similar transition states and bifurcations. 

The flux of ergode through a dividing surface as a function of energy and the effect of the 
Morse bifurcations was considered in [MS 14] and commented upon in Section 3 for planar atom- 
diatom reactions. Seeing as the systems representing these examples have more than 2 degrees 
of freedom, apart from the unrealistic planar atom-frozen diatom case, the Morse bifurcations 
do not have a significant effect on the flux, which varies C m_2 smoothly through these. 

The actual use of capture rates as bounds on reaction rates is debatable, but largely de¬ 
pends on the reaction being considered. However, these were only chosen to provide relatively 
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simple examples of Morse bifurcations, and to show the importance of the attitude and angular 
momentum degrees of freedom in binrolecular reactions. 

We considered binrolecular reactions with Euclidean symmetry and reduced them accord¬ 
ingly. Even though (symplectic) reduction theory is a highly developed subject, we faced a 
number of difficulties when considering these examples. Setting aside the fact that singular re¬ 
duction was required, due to the nature of the rotational symmetry, and that singular cotangent 
bundle reduction is still not a complete theory, there is a large gap between the reduction theory 
literature and applications. Even if one restricts one’s attention to the principal non-singular 
stratum, it is not an easy task to find suitable charts. Some of the literature avoids charts 
altogether, focusing instead on the global geometric properties of the reduced spaces, whereas 
the celestial mechanics literature considers charts for different regions of the reduced space. The 
most common approach in the molecular literature is to restrict one’s attention to non-collinear 
configurations such that the gauge theoretic approach to cotangent bundle reduction provides a 
set of charts, as reviewed in Appendix A. However, here we face the opposite issue, namely the 
reduced space is an S 2 fibre bundle, due to the angular momentum degree of freedom, but the 
global nature of this bundle is generally not discussed in the literature. We feel that more work 
is needed, both on charts for the reduced spaces and on their global nature, and that this would 
improve our understanding of molecular reactions, and also other n-body systems. 

By considering normal molecules, with a fixed equilibrium configuration and energies below 
that at which either of the two molecules dissociates, collinearity becomes a decreasing concern 
with increasing size of the molecules, namely codimension-(2n* —5) where rq > 3 is the number of 
atoms in the ith molecule (i = 1, 2), and the chemistry of the molecule is not taken into account. 
However, for smaller molecules, higher energies, or other transport problems we may need to 
consider collinear configurations. For non-zero angular momentum, collinear configurations are 
a subset of the principal reduced stratum. However due to collinear configurations having non¬ 
trivial configuration space isotropy, we cannot find charts via the gauge theoretic approach to 
cotangent bundle reduction. The issue is therefore not one of reduction per se, but only of 
finding suitable coordinates. The transport problem and bifurcations of transition states will 
be the same as those considered in Section 3. For more than seventy years, chemists have been 
using charts obtained by modifying gauge theoretic cotangent bundle reduction [Say39]. The 
idea is to pass to a rotating frame in which the collinear (equilibrium) configuration is along 
a chosen axis, say the aq-axis, but retain the remaining rotational symmetry (about aq) as an 
internal coordinate. Then by choosing the Eckart convention and the non-gauge invariant form 
of the kinetic energy, we find that the Lagrangian is not a function of the angular velocity about 
the collinear axis uq, so we can obtain a Hamiltonian that is not a function of the first angular 
momentum component Zi. That is, l\ is replaced by the canonical momentum conjugate to the 
“internal” rotation about the aq-axis. These charts were first considered by Sayvetz [Say39], 
though nowadays they are often attributed to Watson [Wat70]. This procedure can be justified 
geometrically by applying the slice theorem (see e.g. [OR04, Section 2.3.14]) to configuration 
space in a neighbourhood of the collinear configurations, and then lifting the charts obtained to 
the cotangent bundle [RSS06]. Actually, with this understanding, charts can be obtained that 
are not those of the Eckart convention, i.e. other gauges and internal coordinates. This was 
used in examples by Kozin et al. [KRTOO]. Note that this is just the splitting of coordinates 
into internal coordinates and rotations and not an actual reduction, cf. Appendix A. With these 
charts, we cannot simply pass to Serret-Andoyer coordinates to reduce the symmetry, seeing as 
the Hamiltonian is not a function of l \. This is generally not addressed in literature. 

The transport problems associated with reaction are usually a lot more complicated than 
the capture ones considered here, so their transition states may undergo a number of different 
bifurcations. One simple example which should display much the same bifurcations as those 
seen here is isomerization reactions involving only one molecule. These will be the topic of a 
future publication. 
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We have concentrated here on how the transition state and dividing surface vary with energy, 
but the exact dependence of the capture transition states on the angular momentum and the 
possible loss of normal hyperbolicity for large values should be considered in detail. Due to the 
high degrees of freedom of these systems, this is not a straightforward task. 

It is also interesting to ask whether the dividing surface method can be extended to consider 
reactions with a varying external held or laser pulse, which need to be modelled as a non- 
autonomous Hamiltonian system; or reactions that are not in the (dilute) gas phase, for which 
the product kinetic approximation leading to a low dimensional Hamiltonian system is not 
valid; or even reactions out of equilibrium. Reviews of the dividing surface approach applied 
to the basic transport scenario of flux over a saddle for non-autonomous systems and Langevin 
systems can be found in [BMH + 08] and [KKllb], respectively. However, more work is required 
to consider general transport scenarios and fully understand transport in these systems. 
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A Charts for reduced n-body systems in non-collinear configu¬ 
rations 

Reduced charts for n-body Hamiltonian systems in the non-collinear configurations region can 
be found by considering the Euclidean action of SE(3) = R 3 tx 50(3) on configuration space. It 
is because one considers the SE( 3) action on configuration space that the configurations must 
be non-collinear in order to avoid coordinate singularities when the action is not free. This 
is a different issue from the stratification of the reduced state space. Note that generally the 
configuration and state space isotropy subgroups for a given Lie group G action are not the same, 
instead we have that G z C G q for z = (q,p). The gauge theoretic approach to cotangent bundle 
reduction is nicely reviewed by Littlejohn and Reinsch [LR97]. They however do not consider 
the final step required to reduce the rotational symmetry and fix the angular momentum. This is 
achieved by introducing Serret-Andoyer 6 coordinates, as explained by Deprit [Dep67] (see also 
[DE93, CW12]). These introduce inevitable coordinate singularities on the angular momentum 
sphere, which is probably why Littlejohn and Reinsch avoid them. We shall briefly review the 
gauge theoretic approach for general n-body systems, and introduce our notation. A specific 
choice of charts for n-body system representing bimolecular reactions is given in Section 3. 

Consider a translation-reduced, rotation invariant n-body systems restricted to the non- 
collinear subset (i.e. the trivial configuration isotropy-type submanifold) of configuration space 
Qid C Q = ]R 3 ( n— P and written in the Lagrangian formalism 

- 71—1 

L(R, R) = -Y^m i \R l \ 2 -U(R), 

1=1 

where R = (R±, ■ ■ ■ , R n -i) are some choice of Jacobi vectors, and rn t the reduced masses. The 
Jacobi vectors are not normalised or mass-weighted. We believe that the mass parameters are 
best dealt with by non-dimensionalising the system. The potential U is assumed to be invariant 
under the action of 50(3). 

Pass from the inertial frame {X\, X 2 , X%} to a convenient rotating frame { x\,X 2 ,x 3 }, which 
will depend on the problem at hand, and write Ri = g (t/>) ■ ri , for i = l,--- ,n — 1, where 

*6 Often also referred to as Andoyer or Deprit coodinates. A nice account of their history is given by Deprit 
and Elipe [DE93]. 
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g £ 50(3) is the rotation parametrised by the Euler angles V> = (V* 1 ,^ 2 , ^ 3 ), and r* are the 
Jacobi vectors in the rotating frame. The rotating Jacobi vectors can be expressed in terms of 
3n — 6 internal coordinates q for Qm/SO{ 3) by specifying ri(q), which is called the gauge in the 
physics literature [LR97]. We are effectively considering a fibre bundle 717 ^ : Qid Qid/ 50(3), 
and q are coordinates for the base space. Then, a(q) = g (ip) ■ ri(q) is a section, and the Euler 
angles -0 = (Vh> i> 2 ,4>3) are coordinates for the fibre, diffeomorphic to 50(3). 

In the new coordinates, the kinetic energy is 

3n—6 3 3n—6 3 

2-Efc — ^ ^ q%Kij (q)qj + 2 EE E 

i,j=1 i,j=1 fc=l ij=l 


where u; is the angular velocity, that is the vector corresponding to the skew-symmetric matrix 
fi(i/>) = g T (ip) 9 (VO) f° r which uj x r = fir, for any 3-vector r. We are therefore considering an 
anholonomic frame (or vielbein ) (< 7 , w) for the tangent space at (q,ip), with uj = ^(ip)ip [LR97, 
Appendix C]. The pseudo-metric K{q) satishes 


n —1 

I<ij(q) = ^m k 
k =1 


dr k (q) 

dqi 


dr k (q) 

d qj 


This is the restriction of the Euclidean metric on the (translation-reduced) configuration space 
Qid to the section a(Qid/SO( 3)), and hence a “pseudo-metric” on the internal space Qi d /SO( 3). 
It is of no importance in gauge theoretic terms, but nonetheless features prominently in the 
molecular literature. The moment of inertia tensor I(q) is given by 


n -1 

J («) = I] m k (r k (q) ■ r k (q)I d - r k (q) <g> r k (q)), 

k =1 


or 


n -1 / 3 

hj{q) = m k y ^r k s(q) 2 5ij - r ki (q)r k j(q ) 

fc=l \s=l 


where (g> is the tensor, or outer, product for which r k <S> r k = r k r^, and the gauge potential A(q) 
associated with the Coriolis effect, which is caused by the coupling term, is 


A(q)=I 1 (q)a(q), 


where a(q) = (ai(q),--- ,a 3n _ 6 (g)) and 


n —1 

a i(q) = x 

k =1 


dr k (q) 

dqi 


Equivalently 


n—1 3 

Aj(g) = IJS (q) e stur k t{q) 

k =l s,t.,u= 1 


dr ku (q) 

dqi 


where I ks (q) are components of I~ 1 (q), and e^ k the Levi-Civita symbols* 7 . 

The kinetic energy is gauge invariant, i.e. independent of the choice of internal coordinates, 
but the individual terms are not (see [LR97, Section IV.A]). It is therefore rewritten, in a 
gauge-invariant form, as 


3 n —6 3 3 n —6 

2 Ek — ^ ^ qiKij(q)qj + EE (uji + A ki (q)q k ) I i:j (q ) (ujj + A kj (q)q k ), 

i,j= 1 *,j=l fc=l 

*7 Recall, the Levi-Civita symbol e;,*, is 1 if ( i,j,k ) is an even permutation of (1,2,3), —1 if it is an odd 
permutation, and 0 if any index is repeated. 
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where the metric K(q) = K(q) — A T (q)I(q)A(q) is an actual (Riemannian) metric on the internal 
space, obtained by projecting the metric on configuration space Qid down to the internal space. 
It is therefore positive definite, but non-Euclidian due to the nature of the space [LR97, Section 
IV.C]. 

Finally, pass to the Hamiltonian formalism. The momenta are found (via the fibre derivative 
of the Lagrangian) to be 

,_ dL{q,ip,q,u) _ t ( w a I \ • \ 

H • / j / j Iij\Q) y^j “ 1 “ AkjyQjQk) ? 

3 =1 k =1 

at / / • \ 3 tl 6 3 

Pi ■— q. — Y + / ,Aij(q)lj, 

* j=i i=i 

where l is the angular momentum in the rotating frame, i.e. I = g T (i/j) ■ L. The Hamiltonian is 
then the Legendre transform of the Lagrangian, namely 

- 3n—6 3 

H(q,i>,p,l) = - Y Y^ Pi ~ A ik{ q )h) K 13 { q ) (jPj ~ A jk{ q )h) 
i,j= 1 k=l 

1 3 

+ 2 Y l i p3 ^ l i + U(q), 

*j=i 

where the potential is a function of the internal coordinates only, due to the assumption of 
rotational invariance, and the Euler angles are ignorable. The symplectic form is 

3n—6 3 1 3 

LO = ^ d% A d Pi+ Y A d lj + - Y Gjkh^juW^kvWd-ipu a dip v . 

i= 1 i,j=l £ 

Alternatively, the Poisson bracket for two smooth functions E, G is 

{E, G} = (d qi Fd Pi G - d Pi Fd qi G) + V* (d^Ffy-G - d tj Fd^G) - e ijk lid h Fd lk G. 

Littlejohn and Reinsch derive this in [LR97, Section IV.D]. 

The momenta p are gauge dependent because of Coriolis term A T (q)l. Passing to gauge- 
independent non-canonical momenta 8 , Vi = pi — Aij(q)lj, simplifies the Hamiltonian to 

3n— 6 3 -. 3 

H(q, v, l) = - Y Y v i Rij ^) v 3 + V(q, l ), V(q, l) = - Y h&Wj + U(q), 

i,j= 1 k= 1 i,J=l 

where V is the effective potential combining the centrifugal term and the potential, and 

3n—6 3n—6 3 

w = Y^ d Qi A d Vi + EE Aij(q)dqi A d lj 
i=l i— 1 j =1 

^ 3n—6 3 

+ 2 S E lj (B ki j(q) T tjuv A ku {<[) Aj /V (q ) ) dqi A dq k 

i,k= l j'=i 

3 1 3 

+ ^ ^jj(-0)d'i/>j A d lj + - Y Gjkh^juW^kvWd-i/iu A dip v , 

i,j= 1 i,j,k,u,v= 1 

*8 Littlejohn and Reinsch call these “covariant shape velocities” and denote them v. We shall use the same 
notation, hoping that it will not lead to any confusion, even though it gives Vi = Kij(q)qj. 
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Figure 7: Transformation to Serret-Andoyer coordinates. {X\, X 2 , X 3 } is the lab frame, 
{xi,X 2 ,xz} the chosen rotating frame and l the angular momentum vector. ©2 = |Z|. 

where we have introduced the Coriolis tensor 

^ijk(Q) — Ay k(Sl) dqjAik(.Q) ^kst Ajs (q) Ajt (q) , 

which is a curvature form on fibre bundle (see [LR97, Section III.G]), and simplifies the equations 
of motion. Effectively, this transformation moves the Coriolis effect from the Hamiltonian to the 
symplectic form, in the second and third terms. It is similar to using non-canonical coordinates 
for a charged particle in a magnetic field, such that the effect of the Lorentz force comes from 
the symplectic form, see e.g. [Mar92, Section 2.10]. The molecular literature usually does not 
pass to the gauge-invariant form of the kinetic energy, see discussion in [LR97, Section IV.F], 

By introducing the rotating frame, we have split the coordinates into internal coordinates q , 
(ignorable) rotations ip and their momenta, but we have not actually reduced the system. Since 
(ip, l ) are non-canonical, the fact that ip are ignorable doesn’t lead to l being constant. We can 
however pass from the non-canonical (ip, l) to canonical Serret-Andoyer coordinates ( 9 , 0) which 
consist of the total angular momentum |Z| plus two projections of l, which we are free to choose, 
and three angles. The choice of projection onto the x± and Ai-axis is shown in Figure 7. 

We immediately note that 


l = *(03,02,03) = (03, y/el-el Sin 03 , ^02 -02(508 03 ), 

which we need to transform the Hamiltonian function, whereas the relations between the new 
angles 0 and the non-canonical angular momentum coordinates is less straightforward and de¬ 
pends on the original choice of Euler angles. These are of no use to us here, but can be found 
in [Dep67] and [DE93], where (0,0) are shown to be canonical coordinates. The Hamiltonian 
in these new coordinates is 


- 3n—6 3 

H(q, 01,02, 03,W,©1,©2,©3) = 2 ^2 V i KV ( ( l) V j + V (<h d 3,02,03)- 

i,j=1 k=1 

Therefore, the system is reduced by eliminating the ignorable degree of freedom (0i,0i), fixing 
02 = A, which is the constant absolute value of the angular momentum, and eliminating 62 . The 
remaining angular momentum coordinates (03, ©3) are the canonical latitude and longitude on 
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the angular momentum sphere § 2 , henceforth denoted z\ = ( q\,p \), and there is a coordinate 
singularity at p\ = A. The reduced Hamiltonian function is 


3n —6 3 


H [q,q\,v,px] A) = i EE ViIC l \q)Vj + V(q,z\; A), 

M=1 fc=l 

1 3 

H(g,^;A) = - l i (z x -,\)ri(q)l j (z x ;\) + lHq), 


*d=i 


and 


3n—6 


3n—6 3 


CU = 


2 dqi Advi+ ^Y1 A ij{q)d Zxk lj(z \; A)d% A d 2 Afc 

2=1 2=1 Jf = 1 

3 n—6 3 


+ 2 ^ ^ ^(z A ; A) (B kij (q) + e JU „H fcu (g)A^(g)) d% A d% + dg A A dp A . 
i,k =1 i=l 

The choice of projection is equivalent to a choice of which axis to use as a longitude for § 2 . 
The transformation for other projections is equivalent. By considering e.g. minor and major 
principal axes, we get two charts that cover the whole of § 2 . 

B Spatial atom-molecule scaling 

This appendix contains the scaling of the moment of inertia tensor, the gauge potential and the 
reduced metric for the spatial atom-diatom molecule capture problem of Section 3. 

The intermolecular distance degree of freedom is scaled as x = £~^ x x and p x = e c p x , whereas 
the intramolecular coordinates are shifted such that b = 0 at the equilibrium, and then scaled 
as b = £ b b and p b = E^Pb- Thus, the scaled, rotating frame Jacobi vectors are 

r n b (q) = £ c l Pn b (x) = £^(1, 0) 

3n b —6 

n{q) = 9b(/3 ) • pi(b ) = g b (P) • (p° + £ft ^ 4j b j) + ° ( £ l ) > * = 1, • • • , - 1, 

3 = 1 

where p® € M 3 are equilibrium configuration vectors, 9b(/3) £ SO(3)/SO(2) determines the ori¬ 
entation of B and the 3(n& — l)(3n& — 6) constants pjj^f,. determine the intramolecular coordinates 
b and shall be chosen along the lines of the Eckart [Eck35] and Sayvetz [Say39] conventions for 
normal and anomalous molecules in the small vibration regime. 

If we write the moment of inertia tensor as I(q) =: I c (q ) + Ip(q), then 

n b -1 n b —i 

I[3 = y, (r k ■ r k I d - r k <g) r k ) = G b { ^ (p k • PkJd ~ Pk® Pk))Gi , 
k= 1 k= 1 

where GbPk = 9 b ■ Pk- Thus, it scales to 

I(q ) = e- c 2 I c {x) + G b (/3)I° b Gl(/3) + O ( e J) 

where I c (x) = £“ 2 mix 2 Diag(0,1,1), and we choose 7° = T)i&g(p b i, Pb 2 , Pb 3 ) with p b i > p b 2 > 
Pb 3 i v i a 9 b(P) and the rotation of B about aq, i.e. g(tjj). The inverse moment of inertia matrix 
exists, since I is a real, symmetric and positive definite for B non-collinear, and scales as 

e° e 2 £ 2 

I ” 1 ~ I £ 2 £ 2 £c | + 

4 4 4 
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For our choice of coordinates and gauge, the gauge potential A(q) has 

a x {q ) = 0, 

n b — 1 rij, — 1 3n b —6 

a b i{q) = G b (P ) Y Pk x Pw + £feG 6 (/3) (plj x /4)6 j = a^(/3) + e b a^(/3,6), 

fc=i fc=i j=l 

a/3i(^) = a%{P)+0{£ b )% 

and we ask that a^(/3) = 0 for all i = 1, • • • , 3n& — 6, i.e. 

n b -l 

Y 0°fc x p ki) = °> Vi = 1, • • • , 3n & - 6. 

fc=i 

This is known as the Eckart condition, and imposes 3(3n.ft — 6) conditions on p l kvj . The gauge 
potential scales to 

/ 0 £° 0 \ 

A (q) ~ I 0 e 2 c 0 + • • • , 

\ o el 0 ) 

where A° 1 (/3) and x,/3) for i = 2,3. 

Finally, for the reduced metric K{q), which is a real, symmetric, positive definite matrix, we 
consider the pseudo-metric K(q ) and write K(q ) = K c (q ) + Kp(q), where 


and 


/ 1 0 0 \ / 0 0 0 \ 

iV c = 0 0 0 , Kp = ( 0 Kp Kg h ] , 

v 0 0 0 / V 0 kj b k b J 


k pij (q) =Y dGbdGb -° 

k =1 
nj,—1 


dp. Pi ‘ ~Qp~ p k + <9 ( £ fe) — Kpij{P) + 0 (eft) 


A>,(g) = 2 ' Gbplj + O (eft) = A^.(/3) + 0 (e 6 ) 

fc=i Pt 
n b -l 

kbij(q) = Y, Pki ' Pkj + ••• = Aftjj + ... 

fc=l 


We ask that = (Up i )~ 1 Sij for all i,j. That is we choose Williamson normal form coordinates 
for the intramolecular degrees of freedom. These are (3??.ft — 5)(3n b — 6)/2 conditions on p klJ . 
Furthermore, we claim that due to Eckart condition K^ bi -{P) = 0 for all i,j, i.e. 

n b — l ny~t 

Y ~d^ p ° k ' Gbp1 ^ = °- 

Let us consider the case with i = 1. The Euler angles fi can be chosen in a number of ways, and 
the rotation matrix written as G b (P) = G\{P\)G 2 {P 2 ), where Gi(Pi ) is a rotation by about 
some axis y t . Recall that the symmetry about x\ has been reduced and G b (P) €E § 2 . Thus 

d Pl G b (0) = d Pl G 1 (p 1 )G 2 (p 2 ) = G l {p 1 )G 1 { 7 k)G 2 {p 2 ), 

where Gi(^) is a rotation about y\ by ^ and simultaneously a contraction in the y\ direction 
to zero. This can be seen by considering planar rotation matrices. Then 

Tib 1 ^b 1 ^b 1 

Y W-WPk • G ^)Pkj = Y Gi{\)G 2 (P 2 )pI • G 2 {p 2 )p\ 3 = Y Gi{\)pI ■ Pkji 

k= 1 Pl k =1 k =1 
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where p\- = G 2 {h)p l kj and 


n \,—1 

Pk x 


k=1 


rib 1 

G 2 (/3 2 ) E pg x = 0, 


fc=i 


by the Eckart condition. Thus K^ bl -{j3) = 0, and the same is true for i = 2. The gauge 
dependent term of K{q ) instead scales to 

/ 0 0 0 \ 

A t IA ~ 0 e° c 0 + • • • , 

\ o o o / 

so 

( l ~ ° 0 \ 

#(«) = 0 K°(P) + F 0 (P) + e 2 c F 2 (x,P) + e 4 c F 4 (x,l3) 0 + O (e b ) 

V o 0 D- 1 ) 

and 

(l 0 0 \ 

K~\q)= 0 K^(P)-s 2 c J 2 (x,/3)+£p A (x,P) 0 + O (e b ), 

V 0 0 D b ) 

by inverting the matrix blockwise, and expanding inverse matrices in formal power series. 
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